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The concept of electric and magnetic field lines is intrinsically non-relativistic. Nonetheless, for 
certain types of fields satisfying certain geometric properties, field lines can be defined covariantly. 
More precisely, two Lorentz- invariant 2D surfaces in spacetime can be defined such that magnetic and 
electric field lines are determined, for any observer, by the intersection of those surfaces with spacelike 
hyperplanes. An instance of this type of field is constituted by the so-called Hopf-Rahada solutions 
of the source- free Maxwell equations, which have been studied because of their interesting topological 
properties, namely, linkage of their field lines. In order to describe both geometric and topological 
properties in a succinct manner, we employ the tools of Geometric Algebra (aka Clifford Algebra) 
and use the Clebsch representation for the vector potential as well as the Euler representation for 
both magnetic and electric fields. This description is easily made covariant, thus allowing us to 
define electric and magnetic field lines covariantly in a compact geometric language. The definitions 
of field lines can be phrased in terms of 2D surfaces in space. We display those surfaces in different 
reference frames, showing how those surfaces change under Lorentz transformations while keeping 
their topological properties. As a byproduct we also obtain relations between optical helicity, optical 
chirality and generalizations thereof, and their conservation laws. 



I. INTRODUCTION 



The concept of field lines in electromagnetics (EM) is 
interesting for various reasons. Historically, Faraday con- 
sidered field lines as the fundamental entities of EM (see, 
e.g., Chapter 3 of Ref. [1 ). Moreover, everyone's first ac- 
quaintance with the concept of magnetic field is, probably, 
through the simple experiment with iron filings that makes 
magnetic field lines around a bar magnet visible. On the 
other hand, by the time one is an advanced student learn- 
ing about Special Relativity (SR) and the role Maxwell's 
equations played in its development, one typically does not 
hear much anymore about field lines. And indeed, the con- 
cept of field lines is highly non-relativistic for two reasons. 
First, electric and magnetic fields transform into each other 
under Lorentz transformations, and second, the points on 
a field line all have the same time coordinate. 

As it turns out, one can define electric and magnetic field 
lines in a relativistically covariant manner only for certain 
types of solutions of Maxwell's equations. This has been 
known for magnetic field lines since a seminal paper by 
Newcomb [2] (and as a demonstration of current interest 
in the covariant description of magnetic field lines, see [? 
]). There he also discusses two related concepts that are 
logically prior to covariance of field lines, namely, identity 
of field lines and motion of field lines. For, in order to 
discuss what a given field line for one observer looks like 
to another observer, one needs to be able to identify the 
field line in the first place. Similarly, even for one fixed 
observer, in order to be able to talk about the motion of a 
given field line one needs to identify which field lines at two 
different times are deemed to be the same. As Newcomb 
derived in Ref. [2] , magnetic field lines can be said to move 
with a (position and time-dependent) velocity v if and only 



if (employing units where c = 1) 

V x (E + vx B) =0. (1) 

He further derived that a covariant description of magnetic 
field lines exists if and only if 

E- B = 0. (2) 

We will indicate how these conditions arise in Sections [XX] 
and[V| respectively. In the special case that Q is satisfied, 
we can set v = E x B / B 2 to satisfy ([!]) identically. In the 
even more special case where E • B = and E 2 = B 2 , we 
have \v\ = 1, i.e., the magnetic field lines move with the 
speed of light. Such fields are called "null fields' and they 
are the subject of study in this paper. 

Not all representations of the EM field are equally well 
suited for discussing and making covariant the concept of 
field lines. For example, the standard representation of the 
magnetic field consists of specifying B(r, t) as a function 
of position r and time t. The magnetic field line at time 
to passing through a given point fo is then a parametrized 
curve /(A) found, in principle, by solving the differential 
equation 

d/(A) = B(f(X),t ) 
dA |B(/(A),t )| 

with "initial" condition /(0) = fb- Instead, we will make 
use of a more convenient representation, which does not 
require solving differential equations. When applied to the 
magnetic field, the representation goes under the name of 
Euler potentials [3] . It leads naturally to a particular rep- 
resentation of the vector potential developed by Clebsch. 
This representation has several advantages in the present 
context: (i) it allows one to directly define magnetic field 
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lines, (ii) it is easily made covariant for certain types of 
fields, (iii) it shows that field lines can be linked only if the 
gradient of one of the potentials is singular, and (iv) in it, 
the magnetic helicity, a quantity known to contain topo- 
logical information about the linkage of magnetic field lines 
[MB], is compactly expressed. These advantages will be ex- 
ploited here as we are interested in displaying topological 
properties possessed by the (magnetic and electric) field 
lines of so-called Hopf-Rahada solutions [7HTT]. These are 
solutions to the free Maxwell equations (without charges 
and currents). In order to achieve the same compactness 
and simplicity in describing their electric field lines, we 
make use of the dual symmetry of the free Maxwell equa- 
tions. In standard notation this symmetry transformation 
is 



E i — y cos BE + sin OB, 
B ^ cos OB -sin OE. 



(4) 



This symmetry, and the geometric conditions a field has 
to satisfy in order to be able to define both electric and 
magnetic field lines, as well as covariance [not to mention 
Maxwell's equations], all turn out to be easily and much 
more compactly expressed within the formalism of Geo- 
metric Algebra (GA) [T2l IT6] . Because this formalism is 
not widely known (yet), we will start with standard non- 
covariant vector calculus and use it to define Euler poten- 
tials and field lines in Section [TTJ and only then will we 
introduce and summarize the most relevant aspects of GA 
in Section [TTT1 That formalism is used to describe EM non- 



covariantly in Section IV (with conservation laws of helicity 
and like quantities easily obtained), followed by the covari- 
ant description in Section [V] All that theory is applied to 
Hopf-Rahada solutions and a description of their field lines 
in Section VI In general, it is not easy to obtain Euler 



potentials for a given solution, but Rahada has provided 
explicit expressions for Euler potentials in his articles [TJ |H] 
about the Hopf-Rahada solutions. In addition, those arti- 
cles contain many deep insights, only some of which are 
exploited here. In Section |YlI| we consider a much broader 
class of solutions to the free Maxwell equations, obtained 
long ago by Bateman [17 , for which a covariant descrip- 
tion of field lines is in principle possible. We give these 
solutions in an elegant covariant form, without, however, 
being able to find the covariant Euler potentials. 



II. EULER POTENTIALS AND THE CLEBSCH 
REPRESENTATION 

A. Euler potentials 

The Euler representation of the magnetic field amounts 
to writing (for historical and mathematical background, 

see [3]) 



B = \/ax V/3, 



(5) 



where a and (3 are two scalar functions of position r and 
time t. (Euler used this representation not for magnetic 



fields, of course, but for incompressible fluid flows, de- 
scribed by a flow velocity Vf with V • Vf = 0.) The ad- 
vantages of representing B by Eq. ^ are that (i) the con- 
straint 



V-£ = 



(6) 



holds automatically if a and (3 are sufficiently smooth func- 
tions of position, (ii) it conforms to the fact that in princi- 
ple two scalar functions should suffice to express the three 
components of B subject to the single scalar constraint 
and (iii) magnetic field lines at a given time t are tangen- 
tial to surfaces of constant a and to surfaces of constant 
/?, and are thus determined by setting a and /3 at time t 
equal to constants, say, 



mt) = a. 



(7) 



Generically, the locus of points f satisfying these two con- 
straints defines a ID curve in 3D space, and that curve is 
a field line of B. 

A disadvantage is that the Euler potentials are far from 
unique. In fact, it is clear that any two functions /(a, /?) 
and g(a,f3) determine the same magnetic field lines (and 
thereby the same B field) as long as the equations 



f(a,p) = /(a ,AO 
g(a,(3) = g(a ,Po), 



(8) 



have unique solutions for a and j3. This will be the case 
whenever the Jacobian for the transformation (a, /?) \-> 

(/,<?) 



j 

is never zero anywhere, 
have 



df/da dg/da 
df/d/3 dg/d/3 



(9) 



More precisely, in that case we 



B = J- X V/ x Vg, 



(10) 



so that new Euler potentials (/, g) are produced by any 
transformation (a, /?) i-^ (/, g) with J = 1. To demon- 
strate the disadvantage of the non- uniqueness, it suffices 
to note that, if only the Euler potentials were unique, 
then field lines could be given an identity straightforwardly 
through the use of (|7|), by identifying a field line at all times 
by the pair of constants (ao?A))- I n turn this would im- 
mediately yield a (position- and time-dependent) velocity 
v of the magnetic field lines. But as it is, one has to do a 
little bit more work to obtain a consistent definition of a 
field line velocity v. Let us assume that field lines can be 
assigned a field line velocity. In particular, assume there 
is a vector field v such that the time evolution of a pair of 
Euler potentials (a, j3) satisfies 



d t a + v • Va = 0, 



(11) 



3 



Then, in order to obtain equations for the physical field B 
we multiply these two equations by the gradients of /3 and 
a, respectively, and then subtract the two. This yields [3 



x B = d t aVI3 - d t l3Va. 



(12) 



We have not encountered the quantity appearing on the 
right-hand side yet, but the well-trained eye will recognize 
that it, together with Va x Vft will form a Lorentz co- 
variant object (see Eq. (78) below). Taking the curl of (12) 
gives us an equation independent of the Euler representa- 
tion: 



V x (vx B) = d t B. 



(13) 



Since we are interested in solutions to the source-free 
Maxwell equations, we can replace the term on the right- 
hand side by —V x E, thus arriving at 



V x (E + v x B) = 0, 



(14) 



which is the (necessary and sufficient) condition Newcomb 
[2] found for a velocity v to be interpretable as a magnetic 
field line velocity. 



B. Clebsch representation and helicity 

Given Euler potentials a and j3 it may seem straightfor- 
ward now to define a vector potential A from which B can 
be derived as 5 = V x i, namely by 



A = -pVa. 



(15) 



There is, however, a subtlety associated with the magnetic 
helicity H of the field, as has been discussed recently in, 
e.g., Ref. [18 , as well as long ago in [19], and by Rahada 
in the context of the Hopf- Rahada solutions in [6 j. The 
magnetic helicity density h m is defined as 



A-B, 



(16) 



in terms of which the magnetic helicity is given as an inte- 
gral over all space 



=lll h - 



(r)df. (17) 

= 0, and hence H m = 0. 
The next question is whether H m is gauge- invariant (h m is 
certainly not!). Under a gauge transformation 5 A = V^, 
we get, assuming we can apply Gauss's theorem, 



The form (15) implies that 



(V^ • B)df 



{^B-e n )dS (18) 



where the latter integral is a surface integral, with e n de- 
noting the unit vector normal to the surface. For fields 
such that \B\ — >• sufficiently fast for r — » oo, with r = |r|, 
we have that SH m = 0. For the solutions we will consider 
the magnetic field does decay sufficiently fast, and so H m 



is gauge- invariant. In particular, it may now seem that, 
in fact, Hm = for all such fields. One interesting aspect 
of the Hopf-Rahada solutions, however, is that we cannot 
so simply apply Gauss's theorem. In particular, one of 
the Euler potentials — and we pick a here — must be multi- 
valued, and its gradient singular. More precisely, a will be 
given in terms of a particular complex function r] (to be 
defined below) as 



1 {Re(v)\ 
a = — arctan - — — — . 

2tt V Im 0?) / 



(19) 



This defines a multi- valued function, whose gradient is sin- 
gular in the locus of points L where r] = 0, since 



Im(77)VRe(77) - Refa) Vlmfa) 
2tt\t]\ 2 ' 



(20) 



Eq. ( 15 ) implies that A would be singular, too, at the locus 



L. However, we can choose a non-singular vector potential 
by adding a gauge-like term that cancels the singularity of 
the first term, viz. 

A = -f3Va + V*, (21) 
provided we pick \I/ multivalued as well, according to 



(22) 



where f3(L) is the value (presumed unique) attained by 
/3 in the locus of points L. If /3 is determined by \r]\ 2 
alone (which is true for the Hopf-Rahada solutions to be 
discussed below) j3 takes on a unique value /3(L) in all 
points of the locus L. The helicity H m can now be nonzero 
by virtue of the multi- valued character of both ^ and a 
(and this singularity is, therefore, necessary in order to 
describe fields with linked field lines [HE]). In fact, we 
have 

H m = JJJ(\7ax V/3) • Wdf 

= /3(L)(/3 max - /3 min )n 2 (a max - a min ) 2 , (23) 

where n is an integer quantifying the type of multi- 
valuedness of a and it counts the number of branches of 
the function a. For the standard Hopf-Rahada solution we 
have n = 1. This expression makes manifest the precise re- 
lation between magnetic helicity and the Euler potentials. 

III. GEOMETRIC ALGEBRA 

Here we will briefly review the subject of Geometric Al- 
gebra (GA), also known as Clifford Algebra. Extensive 
introductions can be found in articles by Hestenes [I2j [13] , 
who has been advocating its use in physics over many 
years, as well as in the two textbooks Refs. [HJ[T5], and in 
a book of lecture notes [16] . 

One could say that the idea of GA is to subsume both the 
dot product of (3D) vectors and their cross product under 
a single vector product that avoids certain shortcomings 
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of the cross product, but it comes with a multitude of ad- 
ditional benefits, some of which will be made use of here. 
The shortcomings of the cross product are that (i) it is 
not associative and (ii) the geometric notion that ax 6 is 
a vector pointing in the unique (up to a sign) direction 
perpendicular to the plane spanned by a and b does not 
generalize to higher dimensions. The vector product of 
GA is associative and its geometric meaning does gener- 
alize to any number of dimensions, for example, to the 4 
dimensions of spacetime. Now it turns out that an elegant 
covariant description is possible using the GA of 3D Eu- 
clidean space E 3 . That is, we do not need the GA of 4D 
Minkowski space, for reasons discussed in great detail in 
the book [14] and articles by Baylis [20-22 (see also be- 
low). We will, therefore, focus here on the GA of E 3 . We 
denote that algebra by C3 (C for Clifford). 

The vector product of two vectors a and b is constructed 
to be associative and to obey distributive laws for addition 
and multiplication. An additional axiom is that the vector 
product of any vector with itself equal the length squared 
of the vector 



aa 



(24) 



The algebra is then built up from the set of all vectors by 
repeatedly taking sums and products of vectors. We can 
in the end distinguish four (always one more than the di- 
mension of the underlying vector space) different geometric 
types of basis elements. First, we have the standard vec- 
tors of E 3 . Second, we encounter scalars, because of (24). 



The third type arises when we consider the product of two 
linearly independent vectors, say, a and b. Their product 
splits into a commuting symmetric scalar term a • b and an 
ant i- commuting anti-symmetric "bivector" term (our third 
type of term) denoted by a A 6, the wedge product: 

ab = ±(ab + ba) + \{ab - ba) = a-b + aAb. (25) 

The wedge product represents an oriented area for the 
plane spanned by the two vectors. Moreover, it gener- 
ates rotations in that plane (both these meanings easily 
generalize to more dimensions). More precisely, if we use 
an orthonormal basis {efi,^} for a given plane, then any 
vector c G E 3 is rotated in the plane by an angle by the 
transformation 



with 



RcR f , 



R = exp(i0e2ei) 



(26) 



(27) 



and where the f operation reverses the order of vectors in 
any product. One should note here that (e^ei) 2 = —1, so 
that we can also write 



R = cos(^O) + sin(|0)e 2 ei. 



(28) 



This shows explicitly that R is a sum of two different types 
of elements, a scalar and a bivector, and such elements of 
mixed type will always be denoted by sans serif symbols. 



A product of three vectors can be expanded in vector 
terms and a "trivector" term, consisting of a product of 
three orthogonal vectors. (This is the fourth and last type 
of term we encounter.) When we take those vectors to be 
unit vectors forming a right-handed frame, then we denote 
the resulting trivector by 

I = eie 2 e3. (29) 

One can easily verify that this entity is the same for any 
right-handed set of orthogonal vectors. Moreover, it satis- 
fies 

I 2 = -1, (30) 

and it commutes with all bivector s and all vectors. It thus 
plays a very similar role as the standard imaginary unit i 
(and it is denoted as such, by z, in many papers and books 
on GA; here we keep a different notation, to make sure we 
remember / is, in fact, a trivector). The geometric meaning 
of / is that it represents an oriented volume, spanned by 
the three basis vectors (this geometric meaning makes it 
obvious there is only one such element in 3D, whereas in 
higher dimensions there exist multiple linearly independent 
trivectors). / is a pseudoscalar as it changes sign under 
parity reversal. Importantly, in 3D, the standard cross 
product is related to the wedge product by 



a Ab = la x b. 



(31) 



In words, the bivector a A b is dual to the cross product 
a x b. More generally, multiplying a vector by I yields a 
bivector; multiplying a bivector with / yields a vector. 

A general element M of the algebra C3 is then a sum of 
the four types. Thanks to the general relation (31 ) we can 
write this sum as 



M = a + bl + a + lb, 



(32) 



where the terms represent the scalar, pseudoscalar, vector 
and bivector parts, respectively, where a, b are scalars, and 
a, b are ordinary vectors. 

The sum of scalar and pseudoscalar parts is denoted by 
(M) s . Similarly, we denote the sum of vector and bivector 
parts by (M) v . A useful identity is 

(MN),= (NM)„ (33) 

for any elements M and N in C3. Note that the analogous 
relation for the vector part does not hold. It is also useful 
to define "real" and "imaginary" parts of elements of C3 
by 

Re(M) = (M + M f )/2; Im(M) = i(M f - M)/2, (34) 

where the f operation, as mentioned above, reverses the or- 
der of vectors in a vector product; equivalent ly, it reverses 
the sign of I in a decomposition like ( 32 ) . 



Another involution, of great use in descriptions of rela- 
tivity, is the Clifford conjugate, which reverses the direc- 
tions of vectors and pseudovectors but leaves the scalar 
part the same. We denote the Clifford conjugate by a bar, 
like so 



M = a + bl - a - lb, 



(35) 



for M given by (32). 
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IV. NON-COVARIANT DESCRIPTION OF 
ELECTROMAGNETICS WITHIN GA 

A. Preliminaries 

For the description of the EM field we will build on the 
ideas in the textbook [14 j, going beyond its treatment in 
various aspects. This Section gives a non-covariant de- 
scription of EM, the next presents the covariant version. 
The results in this Section pertain to general solutions of 
the free Maxwell equations, not just null fields, except in 



the very last subsection IV D 1 



It might be odd to begin this subsection by referring 
to the standard covariant description of EM, but this will 
explain why the non-covariant GA description is so elegant; 
moreover, it demonstrates that the description of EM in 
terms of the Riemann-Silberstein vector [23] (defined as 
E + iB) is elegant for the simple reason that it is in fact the 
GA description with / replaced by i (see Eq. ( [38] ) below). 

The standard covariant description of the EM field is in 
terms of the antisymmetric tensor F^ v (with the indices 
running from to 3). If we first define the following 4 
elements (since they are not all of the same type we use a 
sans serif symbol to denote them) 



then by defining 



eo 1; ei, 2 ,3 ei >2 ,3, 



F 1 = ii^ (e,e u ) v 



(36) 



(37) 



(using the usual Einstein convention of implied summa- 
tion over repeated indices, and using the Clifford conjugate 
(|35|)) we get the central result [14] 



Fx = E + IB. 



(38) 



Although the object on the left-hand side is covariant, 
the split into electric and magnetic fields is observer- 
dependent. In fact, / is observer-dependent, because the 
three spatial unit vectors are. 



B. Source-free fields 

We now limit ourselves to a discussion of source-free 
fields. For such fields we can introduce transverse (and 
thereby gauge-invariant) vector potentials by 



such that 



IC, 



VF . 



(39) 



(40) 



Indeed, when considering separately the four geometric 
types contained in this equation, we find that the scalar 
and pseudoscalar parts give V • A = and V • C = 0, 
respectively. 
E = - 



The sign convention of (39) is such that 
x C . The reason for introducing C is so we 



can define an electric helicity. In analogy to ( 16 ) we define 
the electric helicity density 



h P = C-E, 



(41) 



and the electric helicity H e is then obtained by integrating 
this density over all space. 

We can inductively define a whole hierarchy of fields 
and/or potentials by 



F n +i — VF n , 
for all integer n. For example, 

F 2 = -Vx5 + IVxi 
For free fields they all satisfy 

(d t + V)F n = 0. 



(42) 



(43) 



(44) 



All these equations are then invariant under the duality 
transformation Q, which is expressed simply as 

F o H>exp(-J0)F o . (45) 

This transformation induces the transformation 

F n ^exp(-/#)F n . (46) 

Having defined the dagger operation before as reversing 
the sign of i", we also have 



F^exp(/0)Ft. 



(47) 



This makes it straightforward to construct quantities that 
are invariant under the duality transformation, namely as 
bilinear quantities containing one term F n and one term 
Fj^ and possibly more factors that are field independent. 
Such quantities are discussed in the next subsection. 



C. Bilinear quantities 

Here we give the relations between certain bilinear geo- 
metric algebraic quantities and the more familiar bilinear 
vectorial quantities. Since the elements of C3 do not com- 
mute in general, it is convenient to consider symmetric and 
anti- symmetric bilinear quantities. For example, we have 



F|Fi + FiF| 
f|Fi - FiF{ 

Re(Ff) 
Im(Ff) 



2{E 2 + B 2 ), 
-4E x B, 
E 2 - B 2 , 
2IE ■ B. 



(48) 



We can recognize the energy density of the field, the mo- 
mentum density, and the Lagrangian density in the first 
three lines. The fourth line gives the quantity that should 
vanish if a covariant definition of electric and magnetic field 
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lines is to be possible, according to [2] 
have 



Furthermore, we 



FiFo + FoFi = 2{E>A-B>d) + 

2I(B-A + E-C) (49) 

F1F0-F0F1 = 2I(ExA-BxC) + 

-2(BxA + ExC) (50) 

These two definitions are not invariant under duality trans- 
formations. Such non-invariant quantities are useful to 
distinguish quantities that are conserved because of some 
symmetry from those that are conserved only for partic- 
ular solutions (see below). Moreover, since measurements 
on the free EM field may make use of electric charges, mea- 
surement can break the duality symmetry, and thus it is 
not true that only duality- invariant quantities are physi- 
cally relevant for free fields, as was argued in [24j [25]. A 
duality-invariant version of Eq. (49) is 



Fif; 



i r o 



FjFx = 2(E 



A + B-C) + 
21 (B • A — E • C) 



(51) 



In the second line the total (i.e., magnetic plus electric) he- 
licity density appears as a pseudoscalar. Combining similar 
relations yields the magnetic and electric helicity densities 
separately: 



B-A = (Im(F 1 )Re(F ) 
E-C = (Re(F 1 )Im(F ) 



Re(F )Im(F 1 ))/2, 
Im(F )Re(F 1 ))/2. 



(52) 



We conclude this subsection by noting that duality invari- 
ance and its relation to helicity and the separation of the 
total angular momentum of light into spin and orbital parts 
has become the subject of very recent studies [? ? ]. 



Cs do not commute, in order to calculate the result from 
the gradient acting on the second term, we first move it to 
the front before taking the derivative. That is, the third 
maneuver is to write 



FiFj = -FjF 1 +2(F 1 F+) , 



(55) 



and use this to rewrite 



V(FxF5) 



VFiFj-V^Fi 
VFiFj + F^Fi 



where in the last line we used 

V(Ft ) = -(VF n )t 
All this together yields 

ft(FxFS) = -V-(FxFS 

-M>. 

= -vVfxfS 



(56) 
(57) 



FlFi 



(58) 



This shows that the helicity is conserved for free fields, but 
so is the quantity 



Q = E-A + B-C. 



(59) 



This quantity Q is not that interesting, however, as it is 
just proportional to the time derivative of A 2 + C 2 (which 
we see when substituting the well known relation E = 
—dtA, as well as the less known but similar relation B = 
—dfC), which is itself a conserved quantity. 



D. Conservation laws for (pseudo) scalar quantities 



Helicities and null conditions 



We consider here in some detail how one may show 
within the GA formalism that a (pseudo) scalar quantity 
is conserved. Take helicity as an example. We use three 
handy maneuvers. The first is to explicitly take the scalar 
plus pseudoscalar part of 



dt ( Fi Fj ) = (d t (F 1 Fl) / 

^ VFlF t^_^ Fl( VF )t ( 
(vFiFj^-^F^. (53) 

The second maneuver is to note that 

^(FFj)) =(v(F 1 Fj\ ) , (54) 



because the gradient vector multiplying a scalar term can 
never give rise to a scalar term. The gradient on the lhs 
acts on both Fi and Fj. But because general elements of 



We now wish to consider the magnetic and electric he- 
licities separately. For this purpose it suffices now (given 
the above results) to consider 



dt(fifo) s = ^(FiFo + FoFx)), 
= -2F? - V • (FxFq), , 



(60) 



where above-mentioned maneuvers were used once again. 
This shows that the magnetic and electric helicities are not 
individually conserved. But for solutions for which F^ = 0, 
they are individually conserved. Solutions of Maxwell's 
equations with F^ = are termed "null solutions" for ob- 
vious reasons, and they have been considered in various 
contexts bot h lon g ago and very recently [T71 126H28] (see 
also Section VII). Now taking the time derivative of F 2 
(which is identically zero for a null field, of course) we find 



d t V\ = -2<F 2 Fi) f 



(61) 
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The quantity on the rhs is the difference between two 
higher-order helicities, not related to the field lines of E 
and B, but to those of V x E and V x B. Clearly, it 
vanishes for null fields. (The sum of those higher-order he- 
licities is, in fact, the optical chirality density, which has 
become of great interest recently [24, 25, 29TF33].) 

We can now repeat the story: since the right-hand side 
must be zero, its time derivative vanishes as well. But that 



in terms of the Clifford conjugate (35). The minus sign 



time derivative is now easily found by reference to (60) 



dt(W s = ^(a t (F 2 F 1 + F 1 F 2 )) s 
= _2F|-V-(F 2 F 1 ),. 



(62) 



This way we encounter higher-order versions of the helicity 
differences, (F n+ iF n ) s , as well as higher-order versions of 
the null conditions, = 0, and they are related through 
conservation laws. 



V. CO VARIANT DESCRIPTION OF 
ELECTROMAGNETICS WITHIN GA 

A. The Baylis method 

Special relativity can be compactly described in a co- 
variant manner by using the spacetime algebra [13], i.e., 
the GA constructed from spacetime vectors living in 4D 
Minkowski spacetime. However, perhaps surprisingly, C 3 
can be used just as well for the same purpose. This has 
been demonstrated by Baylis and coauthors in several ar- 
ticles [2QH22] . There is a subtle difference between the 
two descriptions, related to the fact that the spacetime 
algebra has twice as many basis elements as C3 (it is 16- 
dimensional whereas C3 is 8-dimensional) . Physical quanti- 
ties such as the (rest) mass of a particle or the proper time 
of a particle in its own rest frame play two roles: as the 
zeroth component of the energy-momentum four- vector or 
the spacetime position four- vector, respectively, but also as 
the Lorentz-invariant length of the respective four-vectors 
(in any frame, of course). Whereas C3 uses a single repre- 
sentation for either of these two roles, spacetime algebra 
has two different representations, one for each role. Here 
we follow Baylis' lead (see also his textbook [14]) and use 
the simpler algebra C3. 

Any four- vector is represented by a "paravector", i.e., a 
real element of C3. That is, we can write a general par- 
avector as 



P = Po + P, 



(63) 



with po a scalar and p a 3D vector. For example, the 
spacetime four- vector (as defined relative to some arbitrar- 
ily chosen fixed origin of Minkowski space) is represented 
by the paravector 



t + r. 



(64) 



The Minkowski metric is obtained in natural way by defin- 
ing the "square length" of a paravector to be 



here implies that another important example of a paravec- 
tor, the gradient four-vector, takes the form 



d = d t - V. 



(66) 



The scalar product between different paravectors p and q 
is defined as the scalar part of the product pq, i.e. 



(p,q) = (pq) s = i(pq + qp). 



(67) 



A Lorentz transformation preserves, by definition, this 
scalar product of paravectors. It acts on paravectors as 



LpLt, 



where L can be chosen to satisfy 

LL = 1. 

For example, a pure boost is represented as 
L = exp(?Z?/2), 



(68) 



(69) 



(70) 



where w is the so-called "rapidity." Compare this to the 



description of rotations in space, (26) and (27), which are 



of course, also elements in the group of Lorentz transfor- 
mations. Using (69) we also obtain the combined action of 



Lorentz transformations and involutions: 



L f pL, 
LptlA 



(71) 



Given the condition (69) on Lorentz transformations, one 



easily sees that alternating products of paravectors and 
Clifford conjugates of paravectors transform in a simple 
way. When there is an odd number of paravectors and 
conjugates, then the object transforms just as a paravector 
or as the Clifford conjugate of a paravector. When there is 
an even number of paravectors and conjugate terms in the 
alternating product, then it transforms slightly differently. 
For example, an object = pq transforms as 



0^0' = LOL. 



(72) 



We can distinguish two contributions to = pq: 

= (pq) s + (pq)„ , (73) 

where the first term is just the Lorentz-invariant scalar 
product of p and q, and for the second type of term (i.e., 
a sum of vector and pseudovector parts) Baylis uses the 
name "biparavector." 

B. The EM field 

The EM field is described by a biparavector, namely Fi . 
The covariant form of the free Maxwell equations is really 
the same as before: 



PP 



Pi 



P , 



(65) 



<9F! = 0. 



(74) 
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We can introduce a gauge potential by 



(75) 



Now here is an interesting difference between the covari- 
ant and the non-covariant ways of introducing gauge po- 
tentials. Namely, the equivalent of the electric vector po- 
tential C is defined through 

*Fi = {8C) v , (76) 

where C is real, instead of through 

F 1 = (dA) v +l(dZ) v (wrong!) (77) 

which one might have thought would be the covariant form 
of (Bob and p9). 



Since one can pick only 4 linearly independent paravec- 
tors, Fi, like any other biparavector, can be either writ- 
ten as a single product pq — and then it is called a simple 
biparavector-or as a sum of two such products pq + rs. 
In either case, that decomposition is covariant (unlike the 
split in electric and magnetic fields). Given a biparavec- 
tor Fi it is easy to check whether it is simple or not: it is 
simple if and only if Im(F 2 ) = 0. In particular, if a field 
is describable by the covariant generalization of Euler po- 
tentials (for a discussion of covariant Euler potentials for a 
different class of fields, and from a very different perspec- 
tive, see Ref. [34 ) 



Fi = (dad/3) v , 



(78) 



it must be simple. Given ( |76| ) we can introduce covariant 
electric Euler potentials by 



if, = (d\dfi) v . 



(79) 



These two definitions ( |78| ) and ( |79| ) allow us to define field 
lines in a covariant way, as follows. A 2D surface in space- 
time defined by setting a and j3 equal to constants (note 
we use the paravector r = t + r as argument here) 



a(r) a ; /3(r) /3 



(80) 



is an invariant surface, if we stipulate that a and f3 are 
scalar fields, i.e., that they transform under Lorentz trans- 
formations as 



a'(r') = a(r); /?'(r') = /3(r). 



(81) 



(And this we substitute on the left-hand sides of (80).) 
This invariant surface is swept out over time by a given 
magnetic field line. Conversely, a magnetic field line in a 
given inertial frame of reference is then determined by the 
intersection of the invariant 2D surface with a spacelike 
hyperplane, i.e., a hyperplane whose normal is timelike: 
in other words, it is simply determined by fixing the time 
coordinate in the given frame in Eq. ( |80| ). This then gives 
the covariant definition of magnetic field lines, and it agrees 
with that given by Newcomb in [2 . For finding electric field 
lines in a given frame of reference, we go through the same 
procedure starting with ( 79 ) , and end up setting 



A(r) = A ; /i(r) = /i 



(82) 



and fixing the time coordinate in that frame. This proce- 
dure then determines what field lines in one reference frame 
look like in another. Namely, we simply identify magnetic 
and electric field lines by the pair of constants (qq? Po) an d 
(A ,/io), respectively. 

Returning to the definition of helicity density ( [5T] ) from 
the previous Section we see it is not covariant, as it is not 
an alternating product of paravectors and their Clifford 
conjugates. But it can be turned into a sum of two different 
covariant terms by inserting an arbitrary field-independent 
paravector q and its Clifford conjugate (35), like so: 



h(q) 



FiqFj 



FjqFi 



2(F iq Fj 



(83) 



The first term transforms as a paravector, the second as 
the Clifford conjugate of a paravector, and their sum is 
a scalar. We recognize now the helicity density (51) as 



simply proportional to this scalar. For example, by choos- 
ing q = u := (1,0,0,0) it follows that the helicity density 
equals h(u)/2. Given this, it is easy to see that the helicity 
being an integral over all space of the zeroth component 
of four- vectors is invariant under Lorentz transformations 
(just like electric charge is, for instance). In fact, this in- 
variance also follows from the expression Eq. (23), which 



shows the (magnetic) helicity is expressable in terms of the 
scalar Euler potentials. 



VI. NUMERICAL RESULTS 

A. The Hopf-Ranada solutions 

We now apply the results of the preceding Sections to the 
Hopf-Rahada solutions. Those solutions are usually given 
in terms of complex functions 77 and ( through 



B = 



E 



47TZ (1 
1 



1 77 



2\2 



V77 X Vt?*, 

rVCx VC*. 



4ttz(1 + ICI 2 ) 2 

For completeness, we give them here explicitly 

_ Az + t(A - 1) + i(tx - Ay) 
V ~ Ax + ty + i(A(A- 1) -tz)' 
Ax + ty + i(Az + t(A - 1)) 



(84) 



c 



tx-Ay + i(A(A- 1) -tz)' 



(85) 



with A = (|r| 2 — t 2 + l)/2. But by switching variables from 
77,77* to <j) = arctan(Re(77)/Im(77)) and |7?| 2 (and similarly 
for and £*) we get the (non-covariant) Euler form for the 
magnetic field with 



(86) 



and a similar result for the electric field in terms of £. 
Remarkably, the same functions also can be used as the 
covariant Euler potentials. 



The Hopf-Rahada solutions are known to be null fields, 
i.e., they satisfy both E ■ B = and \E\ 2 = \B\ 2 . We 
can, therefore, define magnetic and electric field lines co- 
variantly, and both types of field lines move at the speed 
of light. 

The magnetic helicity of these solutions is obtained 
through Eq. ( 23 ) by noting that 



Pn 



Pn 



1; P(L) = 1; a n 



1, (87) 



so that H m = 1. The same result is obtained in the same 
way for the electric helicity, H e = 1. These results had 
been obtained before in a purely topological way as the 
index of the Hopf map in [7]. 

We consider here two 2D surfaces in 3D space by con- 
sidering the locus of points where a(r, to) = ao at some 
fixed time to, and, similarly, /3(r, to) = Pq. These two sur- 
faces are not invariant. One topological property of these 
surfaces can be immediately obtained from the character 
of field lines: the field lines constitute a tangential vector 
field that is nowhere zero. So it has an index 0, and there- 
fore, by the Poincare-Hopf theorem, can only exist on a 
closed surface if it has Euler characteristic equal to zero. 
That is, the surface, if it is closed, must be topologically 
equivalent to a torus. In Fig. [I] we plot a particular surface 





FIG. 2: Same as Fig. [T] but in a reference frame moving in the 
positive z direction with speed c/2. 







FIG. 3: Same as Fig.^but in a frame moving in the positive 
direction with speed c/2. 



FIG. 1: Surface fio — 0.5 at two different times, t = on the 
left and t = 1 on the right, in the lab frame. Note the axes in 
left and right plots are not quite the same, and the surface has 
in fact moved up, i.e., in the positive z direction. 



P = constant at two different times in the "lab frame," by 
which we simple mean the frame in which the expressions 
(85 ) hold. We pick here and in all similar figures below the 

2 — 1. The topology 



value /3o = 1/2, corresponding to \rj\ 
of the torus is clearly visible. 

Then we plot the same surfaces in two other reference 
frames, both moving at a speed of c/2, in Figs. [2] and [3] 
That is, in those references frames we pick t' = and t f = 1 
(while making the standard assumption that the origins in 
the different frames coincide at t = t' = 0), and we keep the 
constant Po = 1/2 the same. The plots confirm that the 
Lorentz-transformed surfaces stay closed and topologically 



equivalent to a torus. We note that the surfaces are not 
stationary, so the deformations cannot be explained simply 
from just length contraction. 

In Figs. [4}|6] we then plot surfaces a = ao constant, 
again in the same three reference frames as before, and at 
the same times. The constant is picked the same in the 
three figures, namely ao = 0.1/(27r), so that we can say 
we plot the same surfaces as seen by different observers. 
These surfaces are not closed, but the figures do show they 
have one handle. 

In the last set of three figures, [7| [8] and [9] we plot the in- 
tersections of the surfaces shown in the previous six figures. 
That is, we plot the magnetic field lines corresponding to 
the two constants fto — 0.5 and 2-kolo = 0.1. In addition we 
plot in the same figures the electric field line corresponding 
to the intersection of surfaces A = Ao and \i = no with the 
constants chosen the same as for the magnetic case, i.e., 
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FIG. 4: Surface 27rao = 0.1 at two different times, t = on the 
left and t = 1 on the right, in the lab frame. Note the axes in 
left and right plots are not quite the same, and the surface has 
in fact moved up, i.e., in the positive z direction. 



FIG. 6: Same as Fig. [4] but in a frame moving in the positive x 
direction with speed c/2. 







FIG. 5: Same as Fig. ^] but in a reference frame moving in the 
positive z direction with speed c/2. 



FIG. 7: In blue: magnetic field lines in the lab frame obtained 
as intersections of the surfaces plotted in Figs, [I] and [4] The 
other field line (red) is electric, and is obtained with the same 
procedure applied to the surfaces A = Ao = 0.1/ (2tt) and \i — 
/jlo = 0.5. On the left t = 0, on the right t = 1. 



ao and jhq = /3q. 



B. Extension of Hopf- Ranada solutions 

Ranada has given an extension of his solutions. Namely, 
given the complex function 77, we can pick, according to [8 
the following pair of Euler potentials 



1 



/RefrfM 
arctan - — ; — - 
2tt V Im (^ n ) 



(88) 



for any integer n. That is, the function 77 = \rj\ exp(i^) in 
(84) is replaced by the function r] n = \rj \ exp(ni(f)). 



But these solutions are really the same: after all, setting 
r] n equal to a (complex) constant is equivalent to setting 
77 equal to a (different, complex) constant, and so the field 
lines of 77 and 77 n are the same. Ranada does state that 
the magnetic helicity changes by a factor of n 2 , but that 
arises trivially from the simple relations B n = nB\ and 
A n nAi. 

Whereas surfaces a =constant are manifolds (even if 
not simply connected), surfaces a n =constant for n > 1 
are not, and consist of n branches. We plot an example 
in Fig. [lOj where we chose n = 3 at time t = 0, and 
the constant equals 0.3/(27r). Three branches are clearly 
visible. These three branches correspond to three different 
surfaces obtained by setting a equal to three constants 
ao = 0.1/ (27r) and ao±l/3. Note that these branches may 
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2 3 -i 



FIG. 8: Same as Fig. [7] but in a reference frame moving in the 
positive z direction with speed c/2. 



FIG. 10: Plot of surface of constant — 0.3/(27r) at time t — 
in the lab frame. The mult i- valued function a 3 has 3 branches, 
which may intersect at most in points where the gradient is 
singular. The viewing angle is different here than in previous 
figures, so as to make the structure of the corresponding field 
line(s) more clearly visible (see next figure). 






FIG. 11: Magnetic field lines as determined by the intersection 
of the surface plotted in the preceding figure with the surface 
£ = 0.5. 



FIG. 9: Same as Fig. [7| but in a frame moving in the positive ; 
direction with speed c/2. 



intersect at most at the points where \r]\ 2 = 0, the branch 
cut of the function a n (here, the x axis). One can see in 
the plot that at those points the gradient of a n is singular. 
The field line obtained by intersecting the branches with 

That plot shows 



the surface /3 = 0.5 is plotted in Fig. 11 



that this is not a single field line, but three field lines. 



VII. THE BATEMAN SOLUTIONS 

In his book from 1915 [36], Bateman constructs a broad 
class of solutions to the free Maxwell equations that all 
satisfy the null property, = 0. Recall that for such so- 
lutions electric and magnetic helicities are conserved, and 
field lines can be defined covariantly. The description of 
those solutions in the book comes tantalizingly close to the 
description used in the present paper, with one important 



difference: Bateman uses complex Euler potentials, instead 
of real ones. This does negate their advantage in finding 
field lines or expressing the helicity. The complex solu- 
tions satisfy an additional constraint, which may help in 
constructing real Euler potentials, but we have not suc- 
ceeded in doing so. Even though Bateman's solutions are 
not given in covariant form in his book, it is straightfor- 
ward to convert to covariant notation. 

In more detail then, Bateman finds pairs of complex 
functions <pi and 02, such that (using the lower of the 
ambiguous sign he uses in his Eqs. (2) and (10)) 

Fi = (d<f>id<h) v ■ (89) 
Moreover, these functions satisfy the condition 

(d<i>id<j> 2 ) v = 0. (90) 

This constraint implies V\ =0, as can be easily checked 
within the formalism of GA. 

Examples of pairs of functions <pi^ can be obtained by 
choosing <pi = (3b and <p2 = <^#/2, where a# and (3b 
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(conforming to the notation Bateman uses) must be picked 
with the lower sign in his Eqs. (11) or (13) 



b\ = x sin d — y cos $ — t 

h = 7> (x cost? + ys'mi!) + iz), 



(91) 



where d is a parameter that can be chosen arbitrarily, or 



h = r - t, 
x + iy 



2(z + r) 



(92) 



The Hopf-Rahada solutions are, in fact, a special example 
of the Bateman solutions, as has been shown in [26] . Even 
this correspondence, where both Hopf-Rahada and Bate- 
man solutions are explicitly known, is not easy to demon- 
strate explicitly. 



VIII. SUMMARY 

We described a known set of solutions to the free 
Maxwell equations, the Hopf-Rahada linked and knotted 
fields, in a compact geometric way by using Geometric 
Algebra. In particular, we represented the solutions in co- 
variant form as 



Fx = (dad/3) v , 



(93) 



with a and (3 real scalar potentials. This way of writing 
yields magnetic field lines as the intersections of two such 
2D surfaces a = c\q and ft = (3q. We plotted examples of 



these 2D surfaces and their Lorentz-transformed versions, 
confirming that such surfaces must have genus equal to 1. 
Moreover, this description shows that whereas ordinarily 
field lines cannot be linked, they can when a is multi- 
valued and its gradient is singular. The electric field lines 
are obtained by writing, similarly, 



/Fx = (dXdfi) v , 



(94) 



with / the unit pseudoscalar volume element (so that the 
left-hand side corresponds to the dual of Fi). For the Hopf- 
Rahada solutions, we showed explicitly that electric and 
magnetic helicities are individually conserved in time; that 
helicity is invariant under Lorentz transformations, and 
that their chirality vanishes, all as a consequence of their 
"null property" F^ = 0. Finally, we showed that Bateman's 
solutions, which have the null property as well, can be 
written in the form 



Fi = {d<hd<h) u , 

in terms of complex potentials that satisfy 

(d<p 1 d4> 2 ) v = o. 



(95) 



(96) 



The question how these solutions can be rewritten in the 
form (93) remains open. 
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